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Abstract. We present NICMOS NIC3 observations of 
the Galactic globular cluster NGC 288, taken South - East 
at 2.4 times the cluster's half-light radius in the J and 
H bands. We have detected the cluster main sequence 
down to J ~ 25 and H ~ 24. The corresponding luminos- 
ity function covers the range 3 < Mh < 9 and peaks at 
AIh — 6.8. The theoretical tracks of Baraffe et al. (1997) 
at [Fe/H] = — 1.3 give a mass function which is best fit- 
ted by a log-normal distribution with characteristic mass 
of TTic = 0.42 M0 and a standard deviation a — 0.35. This 
result is fully consistent with the global mass function de- 
rived by Paresce & De Marchi (2000) for a sample of 12 
globular clusters with very different dynamical histories, 
thus confirming that near the cluster's half-light radius 
the mass function appears "undistorted" by evaporation 
or tidal interactions with the Galaxy and should then re- 
flect the initial mass function. 
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1. Introduction 

The Galactic globular cluster NGC 288 is located about 
one degree away from the South Galactic Pole, in a re- 
gion of negligible interstellar extinction, E{B — V) = 0.04 
(Alcaino et al. 1997), and very low contamination by field 
stars. The cluster is at a distance modulus of {rn — M)v = 
14.7 and is relatively metal-poor with [Fe/H] = —1.3 (Al- 
caino et al. 1997). Its stellar population is characterized by 
the presence of blue stragglers and by an anomalous hor- 
izontal branch (HB) which is almost entirely composed 
of stars blue-ward of the instability strip (Alcaino et al. 
1997). Such a blue HB may well be induced by the second 
parameter, although the study of Buonanno et al. (1984) 
ruled out that the HB anomaly observed in NGC 288 is 
due to any peculiarity in the cluster's He content, metal- 
licity or age. Buonanno et al. (1984) estimated an age of 
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16 ± 3 Gyr from the Am(^rpo-HB) method, whereas Al- 
caino et al. (1997) have derived 14±2 Gyr from isochrone 
fitting techniques. 

NGC 288 has been observed at infrared wavelengths 
only by Davidge & Harris (1997), who have acquired 3 min 
exposures in J and 3 min exposures in if of a central field 
and 48 min images both in J and K of a region at 140'/ 
West and South of the cluster core. In this way, Davidge 
& Harris have resolved the red super-giant branch of the 
cluster down to the upper main sequence (MS) at iiT < 20 
and < J — K < 1. Their isochrone fitting has been 
performed with the oxygen-enriched (at [Fe/H] — —1.26) 
tracks of Bergbusch & Vandenberg (1992) and has resulted 
into an age in excess of 16 Gyr. 

In this paper, we present the first NICMOS observa- 
tions of a field in the outskirts of NGC 288 taken in the J 
and H bands. The data are described in Section 2 and the 
cluster's colour-magnitude diagramme (CMD) and mass 
function (MP) are presented in Section 3. Discussion and 
conclusions follow in Section 4. 



2. Observations and data reduction 

NGC 288 has been observed with the NIC3 camera of the 
NICMOS instrument on board the HST during the paral- 
lel observations campaign (GO 7811). The observed re- 
gion is located 2.4 times the cluster's half-light radius 
{rfii = 2f25; Djorgovski 1993) away from the centre, or 
l.'6E and 5.'1S. Six images of the same field are avail- 
able, both in the J and H bands (NIC3 filters PllOW 
and P160W, centered at 1.1 ^m and 1.6 ^m, respectively) 
for a total exposure time of 45 min in J and 47 min in H. 
The detailed log of the observations is given in Table 1. 

The images have been reduced using the NICMOS 
standard calibration pipeline: they have been first pro- 
cessed with CALNICA for bias subtraction, dark-count 
correction and flat-fielding. They have then been asso- 
ciated and combined with CALNICB, to remove cosmic 
rays and to increase the signal-to-noise ratio. Photometry 
on the images has been performed with the DAOPHOT 
package. Stars have been identified with DAOFIND, by 
setting the detection threshold at lOu above the back- 
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Table 1. Log of the observations 



Image Filter Exposure time (s) 
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191 


96 
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F160W 


255 


96 
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191 


96 
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F160W 
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ground. We have traced the radial profile of each identified 
object and discarded those which showed a full width at 
half-maximum (FWHM) larger than 2'.' 5, being the typical 
FWHM of a point source 1'.'5 in our frames. This procedure 
has produced a sample of 75 stars. We have also detected 
15 extended objects, whose fundamental parameters (flux, 
position and shape parameters) have been measured with 
S-Extractor. 

The field not being crowded, stellar count-rates have 
been measured in fixed apertures of 5 pixel radius (equiva- 
lent to 1"), and the corresponding background values have 
been determined in a fixed annulus with a radius of 7 pixel 
and a width of 3 pixel. After background subtraction, stel- 
lar count-rates have been corrected for the NIC3 intra- 
pixel sensitivity, using the equations computed by Storrs 
et al. (1999, cf. Table 2) in the case of out-of-focus cam- 
paign data, and for the camera being out of focus, making 
use of TinyTim (Krist & Hook, 1999) which simulates the 
PSF of the NIC3 camera with the precise optics settings 
corresponding to a specific filter and to a specific observa- 
tion date. We have used TinyTim to compute two PSFs 
for each frame, one for our observation date (November 
1997) and one for January, 15, 1998 when NIC3 was in- 
focus (in-focus campaigns were carried out in 1998 Jan- 
uary and June). We have calculated the encircled energy 
for a 5 pixel aperture for each PSF and used the flux ratio 
in-focus to out-of-focus to correct our measured count- 
rates. Finally, we have multiplied the sample count-rates 
by a factor of 1.075 to reconduct them to the values mea- 
sured in a nominal infinite aperture. 

The corrected count-rates (c) have been converted into 
the VEGAMAG photometric system by means of the re- 
lation: 



m = —2.5 log 




where U is the count-rate/flux conversion factor 
(known as inverse sensitivity) and Z is the flux of a zero- 
magnitude star in the VEGAMAG system, provided for all 



NICMOS filters and VEGAMAG bands by the HST Data 
Handbook and by the NICMOS Photomet ry Update (cf. 
http: / /www.stsci.cdu/instrumcnts / nicmog ) . 

The main source of error in our photometry is due to 
the intra-pixel sensitivity effect, that causes up to 30% flux 
variations among individual images. With the correction 
of Storrs et al. (1999) the spread in photometry is reduced 
to 0.1 mag in J and 0.09 mag in H. We, therefore, assume 
an uncertainty of 0.1 mag in both the J and H bands. 

3. Results 

3.1. Photometric completeness 

The dereddened {E{B-V) = 0.04; Alcaino et al. 1997) J 
and H magnitude distributions are shown in Figure 1 for 
the samples of stars and galaxies identified in the observed 
NGC 288 field. Stars populate the - 18 < J 25 and 
^ 18 < H <~ 24 intervals, while galaxies distribute in the 
17 < J 24 and - 19 < 24 ranges. The fact 

that stars and galaxies as faint as J or _ff ~ 24 have been 
detected indicates that the cluster field does not contain 
objects fainter than the 24th mag and, thence, that our 
star sample is complete. Completeness is also secured by 
the very low level of field crowding, which allows full object 
detection by eye. 

3.2. The Colour-Magnitude Diagramme 

The dereddened CMD of the stars in our sample is plotted 
in Figure 2. It extends over the H-hand magnitude interval 
18 - 24 and the 0.1 < J - H < 1.4 in colour. 

We have superposed on the observed distribution the 
expected _ff-band magnitudes and J ~ H colours as com- 
puted for low mass stars at [Fe/H] = —1.3 by Baraffe et 
al. 1997 (solid line). The theoretical track is here scaled 
by the NGC 288 distance modulus (to - M)o = 14.7. The 
IR colours provided by Baraffe et al. have been computed 
adopting the most recent non-grey model atmospheres and 
an improved equation of state for low mass stars. The quite 
good agreement between the observed and the theoreti- 
cal distributions (within an observational uncertainty of 
0.2 mag in colour) confirms that we have indeed detected 
the lower end of the MS of NGC 288. 

Particularly evident in Figure 2 is the change of slope 
occurring at at ~ 20 and J — H = 0.8. This is due 
to the atmospheric opacity being mostly produced by H2 
molecules for stellar masses lower than 0.5 M©. 

3.3. The luminosity and mass functions 

In spite of the limited number of stars in our sample, we 
have been able to trace the luminosity function (LF) for 
the observed cluster field, as shown in Figure 3, where the 
number of stars observed in each 0.5 mag bin is plotted 
as a function of the iJ-band magnitude. The LF displays 
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18 20 22 24 
H 

Fig. 1. Magnitude distributions of the star and galaxy samples detected in the observed field 



the same general features found by Parcscc & Dc Marchi 
(2000) for a sample of twelve Galactic globular clusters: a 
peak at Mh — 6.8 (corresponding to M/ ~ 8.2) followed 
by a rapid decrease towards fainter luminosities. 

Wc have used the mass-luminosity relation correspond- 
ing to the theoretical track plotted in Figure 2 to convert 
the LF into the cluster's MF. Since the observed LF is the 
product between the MF and the derivative of the mass- 
luminosity relation, we have first adopted a specific MF, 
subsequently derived the mass-luminosity of Baraffe et al. 
and computed their product. The latter has finally been 
compared with the observed LF of Figure 3 until a good fit 
was found. As our first attempt, we have assumed a power- 
law mass function in the form of dN / dlog{m) oc (us- 
ing this notation, Salpeter's IMF would have x = 1.35). 
We have found, however, that no single power-law mass 
function can fit the LF of the observed NGC 288 field. In- 
deed, an index x = —0.1 can well reproduce the bright 
portion of the observed LF at Mh < 6.8, but it overesti- 
mates the number of stars fainter than the LF peak (cf. top 
panel of Figure 4). On the other hand, an index x = —1.1 
properly fits the faint section of the LF at Mh > 6.8, but 
it predicts a number of bright stars a factor of ~2 larger 
than observed (cf. bottom panel of Figure 4). 



As an alternative to a power- law, Paresce & De Af archi 
(2000) have shown that a log-normal distribution peaked 
at 0.35 Mq gives a good fit to the MF of globular clus- 
ters for masses smaller than ^ 0.8 Af© (i.e. for stars that 
are still on their MS). Our data show that the LF ob- 
tained for MS stars in NGC 288 is fully compatible with a 
log-normal distribution of the form: 



( ) 


= A- 


log{m/mc) 


\dlog{m)J 




y/2a 



where A is a normalization constant, provided that the 
characteristic mass takes on the value of nic — 0.42 and 
the standard deviation is ct = 0.35 (see Figure 5). 

We note here, as already pointed out by Paresce & Dc 
Marchi (2000), that the fact that a log- normal distribution 
is a viable form for the MF of a cluster for stars with mass 
m < 0.8 Mq docs not imply that the same hmctional form 
is appropriate at higher masses as well, where a power- 
law distribution seems more appropriate (see Elmegreen 
1999). 
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Fig. 2. Dcrcddoncd colour-magnitude diagram. The the- 
oretical track of Baraffe et al. (1997) computed for 
[Fe/H] = — 1.3 has been superposed to the observed dis- 
tribution. Uncertainties amount to 0.1 mag and 0.2 mag 
on H magnitudes and on J — H colours, respectively. 
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Fig. 3. The luminosity function of the field in the H band. 
A distance modulus (m-M)o = 14.7 has been adopted 



4. Discussion and conclusions 

Comparing the best fitting LF of Figure 5 with the results 
published by Paresce & De Marchi (2000) indicates that 
the MF of NGC 288 is fully compatible with the global 
function found by these authors for a sample of twelve 
Galactic clusters spanning a wide range of metallicity, dif- 
ferent distances from the Galactic centre and plane, and 
structural parameters (core and half-light radius and con- 
centration ratio). 




Fig. 4. Two power-law mass functions, dN / dlog{m) oc 
have been fitted to reproduce the luminosity func- 
tion observed for NGC 288. Top panel: an index x = —0.1 
well matches the bright tail of the LF at Mh < 6.8 but 
fails to reproduce the observed number of fainter stars. 
Bottom panel: an index x = —1.1 fits well the faint tail of 
the LF at AIh > 6.8, yet it overestimates the number of 
bright stars by nearly a factor of 2 



NGC 288 is characterized by a highly elliptical orbit 
which makes disruption by tidal shocking quite effective 
and more important than disruption by internal two-body 
relaxation, according to the calculations of Dinescu et al. 
(1999). In particular, Gnedin & Ostriker (1997) have de- 
rived a disruption time of ~ 1 Gyr for NGC 288, which 
suggests that this cluster might have experienced a strong 
interaction with the Galactic tidal field during its lifetime, 
such that it could be totally dissolved within the next Gyr 
or so. 

Using the log- normal MF determined so far, we have 
computed the index A log N defined as the logarithmic 
ratio between the number of stars with mass m = rric = 
0.42 M© and the number of stars with m = O.7M0. 

Paresce & De Marchi (2000) have defined such a param- 
eter in order to quantitatively describe the shape of the 
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Fig. 5. The best fit to the observed himinosity function is 
given by the log- normal distribution ln(dN/dlog(m)) = A 
- {\og{m/mc)/V2a]'^ where nic = 0.42 M© and a = 0.35 



MF and correlate it with the properties of the cluster (see 
their paper for more details). In Figure 6, wc have plotted 
the value of A log N for NGC 288 along with those of the 
clusters studied by Paresce & De Marchi as a function of 
their expected time to disruption. 

We note here that Paresce & De Marchi have measured 
A log N on the global MF of the cluster, i.e. the MF after 
correction of the effects of internal dynamical evolution 
(mass segregation) . The global MF has been shown by De 
Marchi, Paresce & Pulone (2000) to closely approach the 
local MF, provided the latter is measured near the clus- 
ter's half-light radius (rhi)- Our data, however, have been 
taken at ~ 2Arhi, so that a correction would be required 
to our value of A log N for NGC 288. Unfortimatcly, wc 
do not have enough data (surface brightness and radial 
velocity profiles and more than one LF at different radial 
locations) to investigate the dynamical structure of the 
cluster and correct for the effects of mass segregation. On 
the other hand, as shown by De Marchi et al. (2000), at 
^ 2Arhi the MF is expected to be steeper than the global 
MF in the mass range of interest here (0.4— 0.8 M©) and, 
as such, the value of A log N is to be considered here an 
upper limit. 

Thus, the position of NGC 288 at the opposite corner 
of NGC 5272 confirms the scattered nature of the distri- 
bution in Figure 6, hence the absence of any correlation 
between dynamical evolution and the MF observed near 
the clusters half-light radius as pointed out by Paresce 
& De Marchi. Indeed, compared to the bulk of globu- 
lar clusters with a mean disruption time ~ 32 Gyr, 
NGC 288 appears to have suffered the highest degree of 
erosion from the Galaxy. Yet, its MF is nearly identical 
to what observed for the cluster sample of Paresce & De 
Marchi (2000). So, although it has not been determined 
precisely at the rhi , it is most likely that the observed MF 
closely reflect the cluster's initial mass function (IMF). 
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Fig. 6. The plane (Disruption Time) vs. 
log[MF(mc)/MF(O.7M0)] from Paresce & De Marchi 
(2000) where NGC 288 has been added. The position of 
NGC 288 (upper limit in A log N) confirms the scattered 
nature of the distribution so that no correlation between 
dynamical evolution and mass function may be drawn 



An example of the opposite case is offered by NGC 6712, 
a cluster with a dynamical history not to different from 
that of NGC 288, yet with a very different MF. Its MF, 
measured by De Marchi et al. (1999) at 1.7 rhi peaks at 
0.75 Mq and slowly drops down to 0.3 Mq, while the aver- 
age MF of all the clusters in Figure 6 is always increasing 
in the mass range 0.75 — 0.3 M0. The lack of low-mass 
stars (less massive than 0.75 Mq) in NGC 6712 has been 
attributed by De Marchi et al. (1999) and by Takahashi & 
Portegies Zwart (2000) to the effect of stripping by tidal 
interactions with the Galactic bulge. 

According to Pryor & Meylan (1993), NGC 288 and 
NGC 6712 have a concentration ratio of 0.96 and 0.90, 
respectively, so that they have so far experienced a sim- 
ilar degree of internal two-body relaxation. NGC 6712 is 
twice more massive but a factor of 4 closer to the Galac- 
tic centre. The cluster orbits have the same ellipticity 
(0.75) but different orbital periods, namely 230 x 10^ 
yr and 130 x 10^ yr for NGC 288 and NGC 6712, re- 
spectively. Given an apogalactic distance Ra of 11 kpc 
{Rp = 1.8 kpc) and an observed Galactocentric distance 
Rgc = 11.1 kpc (Dinescu et al. 1999), NGC 288 is now 
at its apogalactic point. On the contrary, NGC 6712, with 
Ra = 6 kpc (Rp = 0.9 kpc) and Rgc = 3.5 Kpc, is only 
half-way between its perigalactic and apogalactic points, 
so that it has crossed the Galactic bulge more recently 
than NGC 288. Based on the difference in the orbital phase 
alone, one might speculate that NGC 288 has been able 
to thermalize its mass distribution after the bulge shock 
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while NGC6712 is still suffering from it. This hypothesis 
is, however, not applicable as Gnedin & Ostriker (1997) 
give a relaxation time at half-light radius of 1.4 Gyr and 
0.7 Gry for NGC 288 and NGC 6712, respectively, which is 
a factor of 5 to 6 larger than the clusters orbital period. 
Hence, both clusters have still to reach internal relaxation 
after the last bulge shocking. 

Without more data on the LF of these clusters much 
further out in their periphery, it is not possible to fully 
understand the differences between their MF. It is, how- 
ever, plausible that NGC 6712, given its orbital parame- 
ters, plunges more deeply than NGC 288 into the bulge 
so that the low-mass stars depletion observed in its LF is 
essentially due to stripping by bulge crossing or that the 
severe modification to its MF has been impressed by en- 
counters with molecular clouds existing in the spiral arms 
of the Milky Way where this cluster currently sits. 
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